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Abstract 



Sr2Cu03 and Ca2Cu03 are considered to be model systems of strongly 
anisotropic, spin-1/2 Heisenberg antiferromagnets. We report on the basis of 
a band-structure analysis within the local density approximation and on the 
basis of available experimental data a careful analysis of model parameters for 
extended Hubbard and Heisenberg models. Both insulating compounds show 
half-filled nearly one-dimensional antibonding bands within the LDA. That 
indicates the importance of strong on-site correlation effects. The bonding 
bands of Ca2Cu03 are shifted downwards by 0.7 eV compared with Sr2Cu03, 
pointing to different Madelung fields and different on-site energies within the 
standard pd- model. Both compounds differ also significantly in the magnitude 
of the inter-chain dispersion along the crystallographical a-direction: ~ 100 
meV and 250 meV, respectively. Using the band-structure and experimental 
data we parameterize a one-band extended Hubbard model for both materials 
which can be further mapped onto an anisotropic Heisenberg model. From 
the inter-chain dispersion we estimate a corresponding inter-chain exchange 
constant J± ^ 0.8 and 3.6 meV for Sr2Cu03 and Ca2Cu03, respectively. 
Comparing several approaches to anisotropic Heisenberg problems, namely 
the random phase spin wave approximation and modern versions of coupled 
quantum spin chains approaches, we observe the advantage of the latter in 
the reproduction of reasonable values for the Neel temperature T/v and the 
magnetization niQ at zero temperature. Our estimate of J± gives the right or- 
der of magnitude and the correct tendency going from Sr2Cu03 to Ca2Cu03. 
In a comparative study we also include CuGeOs. 
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I. INTRODUCTION 



Initiated by the discovery of high- Tc superconductivity in cuprate com- 
pounds, there is a renewed and growing interest in the electronic and magnetic 
properties of quasi one-dimensional (ID) structures near half band filling. 
Anionic quasi ID CuOs chains of oxygen corner sharing CUO4 plaquetts are 
present in Sr2Cu03 and in Ca2Cu03 and Cu02 chains of edge shar- 

ing plaquetts are present in CuGeOs Moreover, the first two of these 
compounds can be regarded as 'parent' compounds of more complex struc- 
tures as the double chain compound SrCu02 (which, in notation Sr2Cu204, 
is obtained from Sr2Cu03 by replacing the CuOa chain with a CU2O4 dou- 
ble chain), and of a whole family of celebrated multi-leg ladder compounds 
Sr2nCu2n+204n+2 [@]. A detailed theoretical description of the 'parent' com- 
pounds is naturally a prerequisite for understanding all those materials. It is 
also desirable in view of the related two-dimensional (2D) cuprate structures 
of the high-Tc materials 11 1 , and particularly in view of possible dynamical 
stripe superstructures in the latter [12|. 

Experimentally, Sr2Cu03 and Ca2Cu03 are found to be the best realiza- 
tions of the quasi ID spin-1/2 antiferromagnetic Heisenberg model (AHM). 
Their Neel temperatures, T/v ~ 5K for Sr2Cu03 and T/v ~ 9K for Ca2Cu03, 
are very low compared to the intra-chain exchange integrals Jy ~ 0.2 eV, and 
the ordered moments (<0.1 /iBohr) are extremely small [13-15]. The value Jy 
= 190 meV for Sr2Cu03 jl^ appears to be the record value of an exchange 
integral among all known quasi ID antiferromagnets. The correct description 
of the physics of a magnetic quasi ID system with a weak magnetic inter- 
chain interaction has recently attracted much theoretical attention, see, e.g., 
Refs. 0-11]. For all these reasons, Sr2Cu03 has been announced to become 
a 'superstar' in the field of low-dimensional magnetism in near future pO| . 

In the work of Ami et al. with the assumption of only a very weak 
dipolar inter-chain interaction constant J± ~ 0.01 meV and on the basis of 
RPA spin wave theory, for Sr2Cu03 a Neel temperature below 0.03K [^] was 
conjectured. At that time, due to resolution problems and because of the very 
small ordered moment <0.1/iBohr; neutron diffraction on powders failed to 
detect antiferromagnetic ordering down to 1.5K. However, the relatively small 
inter-chain distances of 3.3 to 3.5A suggest that direct inter-chain hopping, 
which leads to a much stronger kinematic exchange interaction, cannot be 
neglected. The discussion of consequences for the magnetic properties forms 
a main issue of the present paper. 

While numerous band-structure calculations [22| and tight-binding pa- 
rameterizations of one- and multi-band Hubbard model Hamiltonians |23| 
for the quasi 2D cuprate structures can be found in the literature, we are 
aware of only two band-structure calculations for a quasi ID cuprate structure, 
both for CuGeOs |24,25|. On the basis of a wealth of available experimental 
data for CuGeOs large efforts are currently directed to the parametrisation of 
phenomenological antiferromagnetic spin-1/2 Heisenberg models and exten- 
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sions including frustration in the next-nearest neighbor intra-chain exchange 
|26-2£]. To our knowledge, no parameterization on the level of Hubbard-type 
models has been undertaken as yet, although estimates within the Anderson 
impurity model |30| were pointing to strong correlation. The understanding 
of CuOs chain substances is less developed, especially with respect to inter- 
chain interactions, and we shall present here a comparative analysis of both 
cases. 

In contrast to the antiferromagnetic ordering of Sr2Cu03, the less 
anisotropic compound CuGeOs exhibits a spin gap state below T5P=14.2K 
which is accompanied by the occurrence of a period-2 superstructure with a 
very small dimerization amplitude of uq ~ 0.007A. Therefore it has been in- 
terpreted as a spin-Peierls (SP) state. The SP state is supported by frustrated 
second neighbor exchange When doped with Zn for Cu or with Si for Ge, 
also a coexisting Neel state has been found below 4.5K pl| , |32[ . The thereby 
observed magnetic moment, for example of 0.23/iBohr in Cui_2^Zn2,Ge03, is 
significantly larger than the corresponding value of 0.06//Bohr observed in 
Sr2Cu03 suggesting a much larger anisotropy of the latter compound. 

Within the frame of strongly anisotropic three-dimensional Heisenberg 
models, all considered compounds should be described by a dominating intra- 
chain exchange coupling, a small exchange coupling in direction of the shortest 
inter-chain spacing which reduces the strong quantum fluctuations and pro- 
vides a non-zero ordered magnetic moment in the ground-state, and a gener- 
ally very small inter-chain coupling in the third direction to ensure a non-zero 
Neel temperature in accord with the Mermin- Wagner theorem. By apply- 
ing the results of available theoretical approaches to the anisotropic two- and 
three-dimensional Heisenberg models [16|-18,33| one may extract phenomeno- 
logical estimates of the exchange parameters for the chain cuprates, which 
then can be compared to results of electronic structure theory. Interestingly, 
different approaches yield significantly different predictions. The accurate 
determination of the Neel temperature of a strongly anisotropic Heisenberg 
magnet is still an unsolved and challenging theoretical task. We review here 
two frequently used approaches and compare their predictions from our esti- 
mates of exchange parameters with experimental data. 

In the following Sec. II we present band-structure results for Sr2Cu03 and 
Ca2Cu03 which were obtained by applying a linear-combination-of-atomic- 
orbitals scheme to self-consistently solve the Kohn-Sham equation with the 
local-density approximation to the exchange and correlation potential (LDA- 
LCAO). This approach provides us in a most natural way with tight-binding 
parameters. Using these results and experimental estimates for short-range 
correlation effects, in Sec. HI a single CUO3 chain is represented by an ex- 
tended Hubbard model and compared to the situation with CuGe03. The 
magnitude of inter-chain exchange is estimated on this basis. In Sec. IV, 
estimates of Neel temperatures and ordered magnetic moments are derived 
by applying both standard RPA spin wave theory and modern quantum spin 
chain theory. The results are summarized in Sec. V, and an outlook is given. 
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II. BAND-STRUCTURE AND INTER-CHAIN TRANSFER 



The crystal structure of the isostructural compounds Sr2Cu03 and 
Ca2Cu03 is depicted in Fig. 1. Chains of oxygen corner sharing CUO4 pla- 
quettes run along the b-direction. The in-chain Cu-0 bond-length is 1.96 A 
and practically the same in both compounds while the Cu-0 bond-lengths in 
c-direction differ: 1.95 A in Sr2Cu03 and 1.89 A in Ca2Cu03. The two in- 
equivalent oxygen sites are referred to as chain oxygen and side oxygen in the 
following. The shortest inter-chain distance occurs in a-direction and differs 
substantially for both cases: 3.49 A in Sr2Cu03 and 3.28 A in Ca2Cu03. 

The self-consistent LDA-LCAO method has been applied to both com- 
pounds with a minimum basis treating the Cu-(4s, 4p, 3(i), 0-(2s,2p), Sr- 
(5s,5p, 4(i) and Ca-(4s, 4p, 3(i) orbitals as local valence basis states and the 
lower orbitals as core states. The crystal potential is calculated from overlap- 
ping spherical site densities. All basis states are calculated in the spherical site 
contribution to the crystal potential and recalculated in each iteration step. 
The valence basis orbitals have been compressed by an additional attractive 
potential to reduce the overlap among them [^. Due to the relatively open 
crystal structure two empty spheres per unit cell have been introduced with 
empty sphere s and p orbitals at each site. For the exchange and correlation 
potential the parameterization of von Barth and Hedin was chosen and it 
has been calculated in atomic sphere approximation. We show in Fig. 2 the 
band-structure and in Fig. 3 the density of states (DOS) of Sr2Cu03. The 
corresponding results for Ca2Cu03 are similar. The quantitative differences 
between both compounds are discussed below. To check the LCAO-band- 
structure by another method, we performed also calculations using the linear 
muffin-tin orbital (LMTO) approximation. We found no substantial differ- 
ences, only the overall bandwidth of the whole pd band complex was found 
to be somewhat smaller in the LMTO results (see numbers below). 

As expected from simple chemical considerations of covalency, there is a 
single, well separated, nearly one dimensional, half-filled antibonding band 
crossing the Fermi level with large dispersion in b-direction (see Fig. 2). The 
width of this band is about 2.2 eV (LCAO) or 2.0 eV (LMTO) for both 
compounds. The characteristic quasi-lD van Hove singularities near the band 
edges are clearly seen in the DOS (see Fig. 3). A tight-binding analysis of the 
orbitals involved shows that in first approximation the Cu 3d^2_y2 as well as 
the side oxygen 2pz and the chain oxygen 2py orbitals are of direct relevance. 
Only a negligible admixture of Cu 4s states can be detected near both edges of 
this antibonding half-filled band. Its weight as determined by the ratio of the 
corresponding areas under the DOS curves (see insert of Fig. 3) is relatively 
small (0.3 per cent for Sr2Cu03), but it increases to a weight of 2 per cent 
for the Ca compound. 

The metallic behaviour of the LDA band-structure is in sharp contrast to 
the experimental observation of large optical gaps ~ 2 eV which are compara- 
ble to the bandwidth obtained above. This points to the necessity of dealing 
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explicitly with the strong on-site Coulomb repulsion at the copper-site. The 
experimental gap cannot be explained by a spin density wave since it is large 
and persists also above the Neel temperature T/v- Instead we have to antici- 
pate the situation of a charge transfer gap between valence states of mostly 
oxygen character and a copper upper Hubbard band above the Fermi level. 

In analogy to cuprates with Cu02 planes the construction of a multi-band, 
Hubbard-like model Hamiltonian would therefore be desirable. However, it is 
well known that such a Hamiltonian can be projected to an effective one-band 
picture which properly describes the low-energy physics |3^,36|. The existence 



of a well isolated, one dimensional band in the present situation (shown in 
Fig. 4 in more detail) suggests such a possibility all the more. We assume 
that the parameters for the one band description can be determined by fitting 
the band of Fig. 4 to a dispersion of the form 

e{k) = —2ti^LDACOs{kyb) — 2t2^LDACOs(2kyb) — 2t±coskxa 

which yields the values listed in the Table []37| . 

In Fig. 4 the dispersion in the a-direction is clearly visible as an energy 
increase with increasing kx by nearly the same amount both at the bottom 
and the top of the band. This dispersion is present through the band and 
gives a value of t± = 25 meV. To be more accurate we determined t± from 
the dispersion at the Fermi level which is shown in the insert of Fig. 4. The 
corresponding dispersion for Ca2Cu03 is significantly larger, by a factor of 
2.5 (2.0 in LMTO). The smaller lattice constant of Ca2Cu03 leads to an 
increase of the inter-chain overlap of the Cu-d and 0-p basis orbitals, but this 
effect alone is too small to explain the strong enhancement. We have checked 
that the transfer in a-direction goes dominantly via the cation Sr and Ca, 
respectively. The two-center Hamilton matrix elements between side oxygen 
and Ca are two times larger than the corresponding ones for Sr. 

Dispersion in the c-direction is found within the level of accuracy of the 
band-structure calculations only (~ 5 meV). This practically missing disper- 
sion in the c-direction indicates also that the inter-chain hopping in a-direction 
takes place horizontally only, with no zigzag contribution in (111) direction 
from side oxygen to side oxygen. 

The comparison of the band structures of the Sr and the Ca compounds 
shows yet another interesting feature: namely, the bonding bands between 
-8.5 eV and -11.5 eV for Sr2Cu03 are shifted downwards by 0.6 to 0.7 eV 
for the Ca compound indicating a significantly different Madelung field. For 
the same reason the side oxygen in the Ca-compound comes out more nega- 
tively charged although the hybridization of its pz-orbital with the Cu 3d^2_y2 
orbital is slightly increased due to the shortened bond-length. This different 
on-site energies may have an influence on the different charge transfer energies 
as discussed in the next section. 

We also analyzed in the same manner (Eq. ||) the linearized augmented 



(1) 



plane waves (LAPW) energy bands for CuGeOs reported by Mattheiss |24| 



CuGeOs differs from the chain cuprates considered above in the structure 
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of the chains. The Cu02 chains of oxygen edge-sharing CUO4 plaquetts of 
CuGeOa result in a more complex highest antibonding band in which two O- 
p states per oxygen site hybridize with the Cu-d orbital. The more complex 
composition of the crystal also manifests in the band-structure. In particular, 
due to the sizeable inter-chain interaction mediated by Ge and due to the 
presence of two chains per unit cell there are two differently filled split-off an- 
tibonding bands. For our qualitative comparison with the above considered 
chains (CuOs chain), we replace them by one half- filled band for the sake of 
simplicity. The tight-binding parameters of CuGeOs (Table) contain a sig- 
nificantly smaller nearest neighbor transfer integral ~ 0.25 eV and an 
anomalously large next nearest neighbor integral |5C] t2,LDA ~ 0.67ti = 0.18 
eV estimated from Fig. 2 of Ref. |24]. The large difference of the transfer inte- 
grals ti^LDA between chains of corner sharing plaquetts and CuGeOs should 
be related to the efficient 180° Cu-O-Cu hopping for the former (a p — d 
bond) compared with the inefficiency of nono" p — d hopping (nearly 135° 
(45°)) for CuGeOs; for further details see Ref. [38|. This special structure 
explains also the relative large next nearest neighbor transfer integrals ^2 in 
CuGeOs due to the involved effective a Pz-Pz hopping along the chain. The 
inter-chain hopping t± can be deduced from the dispersion in the b-direction 
of the LAPW energy bands of CuGeOs and is of the same order as in Sr2Cu03 
and Ca2Cu03. 

Based on both the available experimental data and the band structure 
information obtained here, a semi-microscopic strong correlation model will 
be constructed which then can be mapped approximately onto a spin-1/2 
Hamiltonian to describe the magnetic properties. This is the objective of the 
next section. 



III. MICROSCOPIC DESCRIPTION IN TERMS OF THE EXTENDED 
HUBBARD AND ANISOTROPIC HEISENBERG MODELS 

A. General relations 

Here we parameterize the well-known extended Hubbard model for one 
single chain with hopping terms to first {ti) and second neighbors {t2)'- 

H — ^ ^ ij {^m,,s^'m+i,s + h.C.^ + — y ^ cli^ gCm^sCm,,—s^m,—s 
m;j=l,2;s rn;s 

\ K \ ^ SiSi+i, 

m;j=l,2 i 

where Um = J2s s'^m,s is the density operator and s denotes the spin index. 
We included in (|^ a small, but non-negligible direct ferromagnetic exchange 
which naturally appears if we map a multi-band, Hubbard-like Hamiltonian 
to a one band model Its necessity and its main effects will be discussed 
below. 
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For the low-energy physics, at half-filhng the extended Hubbard model 
dH) can be replaced to leading order in t/U by a spin-1/2 Heisenberg chain. 
It includes also a second neighbor exchange J2 p3] and reads 



H — Ji"^^ SiSi+i + J2 ^ SiSi+2 



jAF 4tf 4*2 



[/ - ' U-V2' 

where the effective exchange integral Ji of the spin-1/2 Heisenberg Hamilto- 
nian of cuprates is reduced from the predominant antiferromagnetic superex- 
change part by the ferromagnetic contribution (Q) 

Jl « - I K I +J^^. 

Notice that within this approach J2 yields a competitional (frustrated) charac- 
ter to the usually dominant short range antiferromagnetic correlations which 
are established by Ji. That term is especially important for CuGeOs. 

The two main parameters U and ti of the effective extended Hubbard 
model are directly related to the optical gap Eg and the exchange integral 
between nearest neighbors Ji which are experimentally accessible. The fol- 
lowing analysis is considerably simplified if the materials of interest are in 
the strongly correlated limit U > At and excitonic effects at zero momentum 
transfer q are not very strong i.e. U ^ tj > Vj; j = 1,2. The parameter 
sets derived below support such a point of view. We take into account the 
effect of the intersite Coulomb interaction Vi by renormalizing the on-site 
correlation in the form U^ff = U — Vi. Then we may use the optical gap 
Eg obtained from the Bethe-Ansatz solution for the pure ID Hubbard model 
given by Ovchinnikov EH], 



„ 00 , 

ml f - idx 

^9 ~ 



Ueff y sinh(27rtix/C/e//) 
[/e//-4ti + 21n2Ji^^ for Ueff^h; J^^ = Atl/Uef f ■ 



In the strong coupling case Eq. (P) can be transformed to the useful relation 

= 0.5Ji^^ (^1 + Eg/ J^^ + 1- 21^2^ . 

It has been assumed that the smaller parameters t2 and | -f^ |, i.e. the hopping 
to second neighbors and the ferromagnetic exchange in (^) , have no substantial 
influence on the charge transfer gap (but t2 enhances the spin gap in the spin- 
Peierls state). 

The presence of a weak second neighbor exchange can be approximately 
described in some cases by an effective renormalized nearest neighbor ex- 
change integral 



J = Ji - r J2 Ji - J2, 

where r = 1 according to Ref. and r = 1.12 according to Ref. |p5|| . 
Recently, Stephan and Penc predicted a strong narrow excitonic peak 
in the density-density response function N{q,uj) of the EHM in the strong 
couphng hmit at the zone boundary q = n/b: 

UJexiT^/b) = U -Vi, 

provided Vi > 0. 



B. Parameter assignment 

In principle, we can determine Ueff and ti from the experimentally mea- 
sured Eg and J values using Eqs. (§-0) which are presented graphically in 
Fig. 5. In the case of Sr2Cu03, very recently also the narrow excitonic peak 
at the zone-boundary (^) and with it U^jf were determined experimentally 

However, to the best of our knowledge, the available experimental in- 
formation on all three systems is incomplete or contradicting each other. 
For instance, for Ca2Cu03 the charge gap determined from the maximum 
of Qme{u>) is £"3 = 2.1 eV but a direct measurement of the J value from 
the magnetic susceptibility does not exist. Interpreting the midinfrared ab- 
sorption as a phonon-assisted two magnon process a value of J = 255 meV 



was reported [47|. For Sr2Cu03 the experimental J values range from 140 to 
260 meV [^,1^,47|. In the following we shall use 190 meV as a representative 
value. According to quite recent data for this system the one-dimensional 
charge transfer gap £'g(Sr)« 1.9 it 0.1 [49| might be somewhat smaller as 



compared to the Ca compound. Strictly speaking, the optical absorption sets 



in already near 1.5 eV 1 41, 42, 49 1. The elucidation of the observed broadening 
of the expected ID van-Hove singularity in terms of the inter-chain interac- 
tion, quantum fluctuations , disorder, and/or excitonic and other many-body 
effects is a difficult problem beyond the scope of the present paper. 

Taking this situation into account, we use the available experimental data 
and also our band-structure results to derive a consistent parameter set of (|2|) 
{ti, t2, Ueff, Vi and | K |) for each of the three substances, separately. Vice 
versa, the demand of internal consistency weights the experimental informa- 
tion. 



1. Sr2Cu03 



Recent electron energy loss spectroscopy data of Neudert et al. |49] allow 
to determine Ueff = 3.15 it 0.1 eV from the maximum of 9me(a;) at the zone 
boundary (|8|). At the same time Eg was measured to be 1.9 ±0.1 eV from 
the data at small momentum. Similarly Eg = 1.92 eV was found from the 
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Raman resonance energy |11| observed for diagonal in-chain (yy) polariza- 
tion only. Our aim is to derive values for the magnetic coupling constants 
from Eqs. with the use of experimental values of U^/f and Eg. Since 

it turns out that the derived J- values depend sensitively on U^jf and Eg, we 
consider two sort of extreme cases. From Eq. (^) we obtain ti = 0.410eV 
with Ueff =3.15 eV and Eg = 1.8 eV (lower bound). That corresponds to 
jf^=213 meV. According to Eqs. (|^,^) that value has to be reduced by the 
frustrated next nearest neighbor exchange J2 of about 12 meV (corresponding 
to t2 = 100 meV from our tight binding fit) and by the ferromagnetic contri- 
bution I K I before it can be compared with the total experimental exchange 
integral J = 190 it 17 meV ||l^. Thus we can estimate a direct ferromagnetic 
exchange of | |« 11 it 17 meV. The slightly smaller K value compared 
with 35 meV for La2Cu04 obtained in Ref. might be attributed to the 
shorter Cu-0 bond-length of 1.89 A for the latter compound. The so derived 
parameter set is listed in the Table. We derived a second parameter set with 
a considerable smaller value of Ji by taking Eg = 1.95 eV and Ueff = 3.25 
eV near the upper bounds of the experimental results. We obtain from (|5|) 
ti = 0.394 eV and correspondingly jf^^ = 190 meV. Such a parameteriza- 
tion is compatible with the total exchange integral 147lg^ meV derived 
from the magnetic susceptibility data. We note that both parameterizations 



are incompatible with the large J values of 246 meV and 261 meV |47| 
derived from midinfrared optical absorption data (the small differences be- 
tween the latter values arise mainly from the adopted phonon frequency of 70 
and 80 meV, respectively, involved in the phonon-assisted absorption process) 
provided there is no sizeable ferromagnetic second neighbor exchange over- 
compensating the ferromagnetic nearest neighbor contribution | K \ and the 
antiferromagnetic next nearest neighbor superexchange J2 = Ato^/Ueff- Any- 
how, the elucidation of the microscopic origin of the apparent discrepancy 
between the magnetic susceptibility and the midinfrared optical absorption 
data analyzed in terms of the simple nearest neighbor spin-1/2 Heisenberg 
model remains a challenging problem. 



2. Ca2Cu03 

The slightly larger charge transfer gap of 2.1 eV suggests also an enhanced 
C/e//-value in comparison with the Sr-compound. That means that it is again 
difficult to find a reasonable parameterization which is compatible with the 
large J value of 254 meV from midinfrared absorption data. Due to the 
lack of experimental information on the magnetic susceptibility we use in the 



following our theoretical estimate of 160 meV for the J value of Ca2Cu03 |36|. 
Adding a ferromagnetic contribution of | X |~ 30 meV (of the same order as 
for La2Cu04) and a frustrating J2 ~ 10 meV we may find jf^ = 200 meV. Of 
course, in the given case this should be considered as a very rough estimate. 
Then, together with Eg = 2.1 eV, we calculate from (|6|) ti = 418 meV. 
According to (||) that corresponds to U(,ff = 3.5 eV, showing the expected 
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enhancement. On the level of the pd- model the reason for the enhanced 
effective on-site interaction should be traced back to a larger Ae^^ « ^e//- 
It seems to be related to a Madelung effect caused by the difference in the 
lattice parameters of the Sr and Ca compound, respectively. This point of view 
is corroborated by our band-structure calculations discussed in the previous 
section: for Ca2Cu03 the distance between the half-filled antibonding band 
and the filled bonding bands is larger by 0.7 eV compared to the Sr case. 



3. Inter-site Coulomb repulsion and comparison with CuGeO^ 



According to the microscopic calculations of Geertsma and Khomski |38| 
the total nearest neighbor exchange integral of CuGeOs Ji(Ge)=11.6 meV 
can be decomposed into an antiferromagnetic contribution of Jj^^ = 30.4 
meV and into a relatively large ferromagnetic one of | K \= 18.8 meV. Let 
us stress again that the effective nearest neighbor transfer integral ti and 
the on-site interaction Ueff are directly related to the antiferromagnetic part 



only. For CuGeOs charge transfer gaps of 3.66 eV |30| and 1.25 eV |43| 



have been reported |50|. Using Eg = 3.7 eV and jf^ = 30.4 meV, the 
main parameters of the extended Hubbard model come out as ti = 0.187 eV 
and Ueff = 4.34 eV. Interestingly, the latter value nearly perfectly coincides 
with the charge transfer energy A=4.2 eV found out from the XPS data 



analyzed within the Anderson impurity model in Ref. |3C | . Within a pd- model 
the significantly enhanced corresponding Aspd value should be attributed to 
the Ge-cations located near the Cu02 chain oxygens. This point of view 
is supported by the following observations. In compounds like Sri4Cu2404i 
where the Cu02 chains under considerations are surrounded by earth alkaline 
cations the corresponding charge gap is reduced to about 2.8 eV [^]. For that 
compound, Aspd as calculated within the ionic point charge model amounts 



to 3.7eV |52|. 



Comparing the data collected in the Table we suggest that Ca2Cu03 
should be somewhat stronger correlated than its Sr counterpart. With- 
out doubt, the most strongly correlated compound among them all is 
CuGeOs having the smallest transfer integral ti(Ge)=0.187 eV and the largest 
C/e//(Ge) = 4.34 eV. The large ratios U^ff/ti obtained in all three cases 
(7.7(Sr), 8.4(Ca), and 23(Ge)) justify a posteriori the use of Eq. (|5|). 

The difference between ti and ti,LDyi may be explained by a renormal- 
ization of the transfer integral ti by the inter-site Coulomb interaction Vi. 
Within the Hartree-Fock approximation, the correction due to Vi leads to a 
renormalized effective hopping integral ti +pVi with the bond order p ~ 2/tt. 
This renormalized hopping integral can be compared with ti^LDA where the 
inter-site Coulomb interaction is already partially taken into account. From 
ti,LDA = ti + pVi and the data of ti and ti^iDA given in the Table we may 
determine Vi = 0.21 eV for SraCuOs, Vi = 0.16 eV for Ca2Cu03, and Vi=0.1 
eV for CuGeOa (here ti = 0.187 eV as estimated above has been adopted). 
Thus, the inter-site Coulomb interactions Vi ~ 0.2 eV of all three compounds 
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are quite close to each other and fulfill the relation Vi <ti U. Notice that these 
numbers for Vi roughly agree with the corresponding 2D-values 0.11 eV or 
0.17 eV given in Refs. |3^j5^, respectively, and the estimate based on the four- 



band model for CuOa chains |54|: Vi « ndnpYpd + n'^{Vdd — V2^dd) ~ 0.23 eV, 
where typical occupation numbers Ud ~ 0.7, Up « 0.13 and Vpd = 1.2eV for the 
copper oxygen inter-site Coulomb interaction and Vi^dd = 0.5 (y2^dd = 0.25) 
eV for the (next) nearest neighbor copper copper inter-site Coulomb interac- 
tion have been taken. 

The value for t2 in the Table was either taken from the fit to the band- 
structure data (= t2^LDA for Sr2Cu03 and Ca2Cu03 ) or inferred from the 



experimentally known value for J2 = 4.3 meV [26| for CuGeOs using Eqs. (5). 



^. Inter- chain exchange 

Now, we would like to give a first estimate of the magnetic couplings 
between chains. The inter-chain exchange interaction J± in the a-direction 
for Sr2Cu03 and Ca2Cu03 (corresponding to the b-direction in CuGe03) will 
be approximated by 

Ji. = (9) 

where we assumed for simplicity the same inter-site Coulomb interaction Vi 
within the chain and perpendicular to it. The corresponding values are listed 
in the Table. The discussion above about a possible direct ferromagnetic 
exchange which leads to a systematic reduction of exchange integrals sug- 
gests that these values should be considered as upper bounds. In the case of 
CuGe03 the so determined J± = 1 meV can be compared with experimen- 
tal data from neutron scattering |55,27| J± ~ 1.1 meV showing a reasonable 
agreement. Notice that in the case of Sr2Cu03 our inter-chain interaction ex- 
ceeds the dipolar interaction evaluated in Ref. by two orders of magnitude. 

The magnitude of the weakest interaction J±^c in the c-direction is difficult 
to estimate theoretically for several reasons. It was already discussed that a 
band-structure hopping integral cannot be given at present. In any case it may 
be expected that J±^c is smaller than the other exchange integrals by several 
orders of magnitude and is difficult to calculate in any case. Treating J^^^ 
therefore as a purely phenomenological parameter in the following, we will use 
to be specific (if necessary) the same value for Ca2Cu03 as has been evaluated 
in Ref. |§ for Sr2Cu03, adopting the dipolar interaction for J±^c ~ 10~^ meV. 

IV. SOME ASPECTS OF THE NEEL STATE 

The magnetic properties of undoped cuprates (i.e. one hole per Cu site 
in the standard pd- model) are usually described by the anisotropic spin- 1/2 
antiferromagnetic Heisenberg model 
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H — 7 ^ JijSiSj , 

<i,j> 

with Jij = J\\{= Ji) for {ij) beeing nearest neighbors in the chain direc- 
tion (that is the b-direction for Sr2Cu03 and Ca2Cu03) and J± for nearest 
neighbor copper sites in a-direction (see Fig. 1). The weakest interaction wiU 
be denoted here by J_l,c- According to the results of the previous sections 
and from the experimental data we know that Sr2Cu03 and Ca2Cu03 are 
characterized by very anisotropic interaction strengths 

^11 > Jl >| Jl,c I, 

the anisotropy being about three orders of magnitude for each inequality. In- 
stead of the spin-Peierls system CuGe03, in this section for comparison we will 
consider the doped compound GeCui_a;Zn2^03 (x=0.034) |^J3^ which shows 
antiferromagnetic order. This is an example for an anisotropic Heisenberg 
problem with weaker anisotropy than Sr2Cu03 and Ca2Cu03. For simplicity, 
we will use for GeCui_^Zna;03 the same parameters which were derived in the 
previous section for CuGeOs. We also neglect here the frustrated exchange 

In the following we review several approaches for such anisotropic systems 
where quantum and thermal fluctuations become important. We will mention 
the usual spin wave approach in self-consistent random phase approximation 
(RPA-SWA) where all directions are treated on an equal and simple footing, 
and the coupled quantum spin chain approach (CQSCA) which involves a so- 
phisticated treatment of the intra-chain direction and a mean-field treatment 
of the remaining perpendicular inter-chain interactions. 



A. RPA spin wave theory 

The RPA-SWA yields simple analytical expressions for the Neel temper- 
ature Tjv and for the magnetization < S\ >= tjiq at zero temperature (see 



[56 1 and references therein). Both quantities can be derived from the 



expression 



2miT) = ^—, V = -yf^cothf^')-l 



1 

where 



il(g) = \JuJq— iofiq) , iOi = 4m(r)J|| {cosQy + Rcosqx + RcCosQz) , 



with Wo = 4m(r) J|| + + J_L,cj , R = J±/J\\ and Rc = J±,c/J\\- We 
put ks = 1- The Neel temperature is defined by the condition of vanishing 
magnetization which yields 

2Tn = JJl{R,Rc), 
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where 



TT^ J J Jo R{1 - COS Qx) + Rc{l - COS Qz) + (1 - COS Qy) 

Expanding I(R, Rc) for Jy ^ J± ^ J±^c gives the approximate expression 

I{R, R,) = ^ (1 + 0.24(ln(i?/i?e) - R)) , 
V R 



which determines the Neel temperature together with (|T^). The zero temper- 
ature magnetization niQ is in the same hmit given by 



0.303 

mo 



1 - 0.386 In (i?) ' 

where the small parameter Rc turned out to be irrelevant. Notice that the 
RPA-description adopted reveals a vanishing magnetic moment in the R 
limit. Thus it differs from the ordinary spin- wave theory which yields a 
diverging expression mo =| 0.5-|-(l/7r)lni? | in the weak inter-chain coupling 
limit. 

Let us now check the above expressions using the estimates of the last 
section and compare them with the experimental data. These data for Tjv and 
the magnetic moment ^^^^ = ghf^Q are given in the Table. In the following 
we will adopt a typical cuprate Lande factor gi k, 2.1 for Cu"^^ p. Using the 
values J|| and J± from the Table and Jj_^c = 10~^ meV, we find Tif = 38 K, 
T^" = 75 K, 11^'' = 0.20HB and /i"^" = 0.26/iB for the Sr and Ca compound, 
respectively. The ratio of the two experimental Neel temperatures agrees 
approximately with the RPA-SWA prediction 



where the logarithmic corrections in Eq. (13) can be neglected since they are 
not very important for the above ratio. However, the absolute values of T/v 
within the RPA-SWA disagree with the experimental data. In the case of the 
more isotropic GeCui-^^ZnajOs we find /i*^*^ = 0.32^B; i-e. a better agreement. 
But also here, the magnetic moment is overestimated by the RPA-SWA. In 
this case that may be ascribed to the effect of the frustrated second neighbor 
exchange. 

For the Sr and Ca compounds one could try the opposite procedure using 
the given experimental data (including Jy) to determine an "empirical" Jj^^. 
Then one finds values for Jj^^ which are two (from T/v) or more than four 
(from mo) orders of magnitude lower than those estimated in the previous 
section. This seems, therefore, to be unrealistic. Despite the fact that it 
gives the correct limits for mo both for R ^ and in the 2D isotropic case for 
— > 1, the RPA-SWA seems to overestimate mo for large anisotropy (i? <C 1) 
quite considerably. That points to the necessity for an improved method. In 
the case of smaller anisotropy (e.g. R ~ 0.1 like for GeCui^j^ZnajOs), the 
RPA-SWA seems to give more reliable results. 



13 



B. Coupled quantum spin chain approach (CQSCA) 

Adopting Schultz's interchain RPA-expression (Eq. (7) of Ref. [|16|), we 
replace 0.5{J± + J±,c) as would be suggested by our strongly "or- 

thorhombic" parameter regime Jy S> J_l S> J±^c- Then one arrives at 

mo = jVR, 

where the proportionality factor 7 is 0.72. (A similar factor 2/7r=0.637, was 
obtained by Fukuyama et al. fS^). The corresponding values for /x'^'^'^ = 
gimo are listed in the Table. 

Analogously, within these theories one expects T/v ~ J± in partic- 
ular, the slightly modified implicit expression for the transition temperature 
proposed by Schulz reads 

r^ = ^JxlnV2 (aJii/Tjv), 

where A ~ 5.8. 

Prom a principal point of view (Mermin- Wagner theorem), it is clear that 
Eq. ( pO[ ) overestimates T/v because it does not depend on J_l,c- However, if 
one assumes that its influence can be described by logarithmic terms like in 
Eq. ( p!6[ ) which then ensure a finite Tat, the relative changes might be quite 
small. 

Like in the RPA-SWA, our estimated values for J± and the experimental 
J|| lead to too large values for T/v and mo. But now, using the experimen- 
tal rriQ and Jy we can determine from Eq. (16) an "empirical" J^*"^ of the 
CuOa chain compounds which is of a similar order of magnitude to our esti- 
mates. The value of is smaller by a factor between 2 and 3 (Sr2Cu03, 
GeCui_2,Zna;03), or 6 (Ca2Cu03) compared to the theoretically estimated 
values (see Table). The Neel temperature can also be used to determine 
which gives similar values showing the internal consistency of the CQSCA. 
But one should keep in mind that Eq. (18) does not fulfill the Mermin- Wagner 
theorem. In that respect we note here an alternative approach to the strongly 
anisotropic Heisenberg model (Castro Neto and Hohn [^) in which the Neel 
temperature was found to depend linearly on J±^c- Naturally, the elucidation 
of the correct description as to how this smallest interaction parameter does 
affect the finite temperature properties remains a challenging problem. With- 
out its generally accepted solution it makes no sense to discuss the absolute 
values of the Neel-temperature beyond an order of magnitude accuracy. 

One possible explanation for the reduction of J'^p in comparison with 
our estimated J± could be the proximity of a spin-Peierls state. Phase fluctu- 
ation effects beyond the mean-field inter-chain approach used in deriving Eq. 
( |l^ ) can then become quite important. Following the renormalization group 
approach of Wang ||l^ for a plane of weakly interacting chains at T = 0, one 
finds a strongly renormalized magnetization which can be traced back to a 
renormalized exchange integral. If that is true, one should expect Ca2Cu03 
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to be much closer to the SP phase transition point than Sr2Cu03. Further- 
more for small exchange integrals compared with the phonon frequency (~ 10 
to 20 meV), phonon exchange gives rise to a quasi-instantaneous interaction 
between localized spins, leading to a renormalization J — > J^jf < J |59|. 

Another possible origin for the difference between J± and might be 
our simple procedure to estimate J± based on the extended Hubbard model. 
It was already mentioned that such a procedure has the tendency to overesti- 
mate the exchange integrals which becomes already apparent for Jy . Last not 
least, one should keep in mind an uncertainty of the band- structure methods 
with respect to transfer integrals as small as in the considered case. The re- 
placement of the full potential in the region in between the chains by empty 
spheres as explained in Sec. II might effect the transverse tails of the Wannier- 
functions which determine the value of the transfer integral t± . A similar effect 
is expected if the exchange and correlational potential is supplemented with 
gradient terms reflecting the strong change in the electron density moving 
from one chain to its neighbors. Anyhow, one would expect roughly the same 
accuracy for the Sr and the Ca compound. In this context the much stronger 
deviation of the magnetic moment of the Ca compound might be related to 
somewhat reduced accuracy of Eq. (16) in less anisotropic cases. 



V. SUMMARY 

Band-structure calculations for Sr2Cu03 and Ca2Cu03 show in addition 
to the expected large dispersion along the chain direction also a remark- 
able dispersion in the crystallographic a-direction, i.e. perpendicular to the 
plane containing the corner-shared CUO4 plaquettes which form the CuOs 
chains. The corresponding inter-chain transfer gives rise to antiferromagnetic 
exchange integrals J± in the meV range. Together with a small dipolar ex- 
change in the third direction J±^c it explains the antiferromagnetic order in 
terms of an anisotropic Heisenberg model. The larger value of J± for Ca2Cu03 
corresponds with the larger Neel temperature and the larger magnetic mo- 
ment in comparison with Sr2Cu03. However, our rough estimation of 
exchange integrals, based on the one band assumption, seems to overestimate 
the difference between both substances. CuGe03 is different from these two 
substances by a much smaller exchange in chain direction and a compara- 
ble large frustration parameter which suppresses the antiferromagnetic state 
and stabilizes the spin Peierls state. The smaller anisotropy becomes appar- 
ent in GeCui_2:Zn^03 where a Neel state was found with significantly larger 
magnetic moments than in Sr2Cu03 or Ca2Cu03. 

The copper oxygen chains of the three compounds under consideration 
can be described within an extended Hubbard model supplemented by ferro- 
magnetic contributions to the nearest neighbor exchange integral. For the Ca 
and Sr compounds excitonic effects in the limit g — > are expected to be weak 
due to small inter-site Coulomb interaction Vi. But near the zone boundary 
q K, Ti/b strong excitonic effects are expected in the framework of the theory 
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developed recently by Stefan and Penc |4^. If our proposed parameterization 
is correct, for Ca2Cu03 (rsS.S eV) and CuGeOs 4.3 eV) these excitonic 
peaks should be observed at higher energies than for Sr2Cu03 (3.15 eV). 

The LDA band structure calculations yield useful insights into important 
material dependent parameters as inter-chain electron transfer and tendencies 
of the crystal field (Madelung ) potential, varying from one substance to the 
other, albeit that the estimate of the on-site and inter-site Coulomb interac- 
tion requires more sophisticated methods such as LDA-calculations with local 
constraint. 

According to our findings, Ca2Cu03 is expected to be a typical ID charge 
transfer insulator analogously to the 2D model system Sr2Cu02Cl2 although 
the inter-chain interaction is intermediate between the strongly correlated 
CuGeOa and the most ID Sr2Cu03. Possibly, the latter system is the weakest 
correlated one of the three. The compounds Sr2Cu03 and Ca2Cu03 offer 
also the opportunity to study in detail the effect of the inter-chain interaction 
provided it can be changed in a controlled way. Indeed, the study of the 
magnetic properties of the alloy system Sr2_xCaa;Cu03 gives an interesting 
possibility to change continuously the magnitude of the inter-chain coupling. 
This is also interesting from the theoretical point of view since it gives a 
possibility to check in more detail sophisticated theories for weakly coupled 
quantum spin chains. 
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TABLES 



TAB. Model parameters for Sr2Cu03, Ca2Cu03 and CuGeOs. The LDA- 
LCAO derived tight binding parameters in the first group of rows are ex- 
plained in Sec. II. The LDA-numbers for CuGeOs are estimated from Figs. 
2 of Ref. 1 24] and Ref. [25|. The second group contains experimental val- 
ues (in the case of several date we prefer the underlined) which were used 
in addition to the band-structure information to estimate the corresponding 
parameters of the extended Hubbard model as well as the exchange integrals 
of the anisotropic Heisenb erg-model (third group of rows). They are derived 
and discussed in Sec. III. The experimental magnetic moment ^u^^^ (given to- 
gether with T/v in group IV) may be compared with ^'-'^^ derived from Eq. 
(19) using the experimental data for the in-chain exchange integrals Jy and 
our estimation of J^. Vice versa, the experimental /U^^^^ determines via the 
same Eq. (19) the empirical inter-chain exchange integrals Jj^^. 
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prrOUD 


quantity 


Sr2Cu03 


Ca2Cu03 


CuGeOa 
(GeCui-^Ziia-Os) 


I 


ti,LDyl/meV 


550 


520 


250 






100 


100 


81 




t±/meV 


20 to 30 


50 to 65 


25 to 33 


II 


Eg/eY 


(1.8 to Lg)^'*" 


(2.1)'= 


(1.25)^1, 




J/meV (= J||/ meV) 


(140)^ (190)^, (260)^ 


(160)^ (254)h 


(11 ± l)j, (22 )^ 


III 


ti/meV 


410 


419 


187 




i2/nieV 


100 


100 


90 






(3.15)^ 


3.5 


4.34 (4.2)*^ 




Fi/eV 


0.21 


0.16 


0.1 




1 K l/meV 


11 


30 


19 




J_L / meV 


0.5 to 1.1 


2.9 to 4.3 


0.6 to 1, (1.1)' 


IV 


Tn/K 


(5)' 


(8.. .10)' 


(4.5)*^ 






(0.06±0.01)' 


(O.OOzbO.Ol)' 


(0.23)-' 


V 


,,CSC /,, 


0.08 to 0.11 


0.19 to 0.24 


0.35 to 0.45 




J77meV 


0.3± 0.1 


0.6±0.1 


0.27 



Raman resonance, Ref. [11| 
^ EELS, Ref. El 



opt. absorption, Ref. H 



^ XPS, Ref. Hi 
XPS, Ref. H 
^ magn. suscept., Ref. ||8|,P] 
s magn. suscept., Refs. [|lO| 
^ midinfrared, Ref. |47| 
' theory, Ref. |6| 
j INS, Raman, Refs. [||j2|j3l|,i5| 



^ Raman, Ref. ^ 
1 ^SR, Ref. m 



INS, Ref. |5| 
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FIGURES 



FIG. 1 Crystal structure of Sr2Cu03 

FIG. 2 LCAO energy-bands near the Fermi level for Sr2Cu03 along high- 
symmetry Brillouin-zone directions within the {kx,ky) plane. The momenta 
are given in units of (tt/o, 7r/b). Strong dispersion can be seen along (0,0)-(0,l) 
(b-direction, parallel to the CuOs chains) whereas a small, but non-negligible 
dispersion in the perpendicular a-direction can be seen. 

FIG. 3 Density of states N{E) for Sr2Cu03. The insert shows the partial Cu 
3d and 4s density of states of the nearly one-dimensional band crossing the 
Fermi level. 

FIG. 4 Dispersion of the nearly one-dimensional band. The insert is for fixed 
ky = kyPgrmi ~ 7r/25. 

FIG. 5 Dependence of the transfer integral ti (upper panel) and the on-site 
Coulomb interaction Ueff (lower panel) of the Hubbard model according to its 
Bethe ansatz-solution vs. inchain superexchange integral jf^ (lower panel) 
in the strong coupling limit for typical values of the optical gap Eg. The ex- 
perimental values for Jf^^ are depicted by arrows. They are determined from 
the total exchange integral J adopting ferromagnetic and second neighbor 
contributions discussed in the text. 
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